Constraints on the density perturbation spectrum from primordial black holes 



OS 
On 

u 

< 
in 



in . 
■ 

o : 
t> . 

On ■ 

^ ' 
Oh; 

6 ■ 

a: 



13 



Anne M. Green and Andrew R. Liddle 
Astronomy Centre, University of Sussex, Falrner, Brighton BN1 9QH, U. K. 

(February 5, 2008) 

We re-examine the constraints on the density perturbation spectrum, including its spectral index 
n, from the production of primordial black holes. The standard cosmology, where the Universe is 
radiation dominated from the end of inflation up until the recent past, was studied by Carr, Gilbert 
and Lidsey; we correct two errors in their derivation and find a significantly stronger constraint 
than they did, n < 1.25 rather than their 1.5. We then consider an alternative cosmology in 
which a second period of inflation, known as thermal inflation and designed to solve additional relic 
over-density problems, occurs at a lower energy scale than the main inflationary period. In that 
case, the constraint weakens to n < 1.3, and thermal inflation also leads to a 'missing mass' range, 
10 18 g M < 10 26 g, in which primordial black holes cannot form. Finally, we discuss the effect 
of allowing for the expected non-gaussianity in the density perturbations predicted by Bullock and 
Primack, which can weaken the constraints further by up to 0.05. 
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I. INTRODUCTION 

Primordial black holes (PBHs) arc formed in the early 
universe if density perturbations are sufficiently large, 
and provide a useful probe of the primordial power spec- 
trum over a wide range of scales. The data from COBE 
and from large-scale structure observations constrain the 
power spectrum on large scales (from about one mcga- 
parsec up to thousands of megaparsecs), whereas PBHs 
may form over a wide range of smaller scales (1 Mpc 
to 1CU 16 Mpc). Limits on their production can be used 
to constrain inflation models in which the perturbations 
grow as one moves to shorter scales, the so-called blue 
spectra. 

There are a number of well-known limits, covering var- 
ious mass ranges, on the maximum allowed mass frac- 
tion of PBHs . Some are imposed at the present 
epoch and some at earlier stages such as nucleosynthesis. 
These constraints fall into two categories, those from the 
effects of Hawking radiation and those from the gravita- 
tional effects alone. The evaporation of PBHs via thermal 
emission has potentially observable astrophysical conse- 
quences, and while no unambiguous detection has been 
made, observations have placed limits on the maximum 
mass fraction of PBHs allowed at evaporation. PBHs 
with mass M < 5 x 10 14 g will have evaporated be- 
fore the present epoch. PBHs more massive than this 
will not have experienced significant evaporation, and 
their present density must not overdose the universe: 
^pbh,o < 1- An additional, less secure, constraint on 
light PBHs can be obtained if one supposes that black 
hole evaporation leaves a stable relic, normally assumed 
to have a mass of order the Planck mass, rather than 
evaporating to nothing. 

It is thought that PBHs can only form in the early 
universe, and that they do so at a time when the horizon 
mass equals the black hole mass. A number of formation 



mechanisms are possible the simplest being forma- 
tion from large density perturbations. The horizon mass 
in a radiation-dominated universe with temperature T is 
given by 



10 18 g 



10 7 GeVS 2 



T 



(1) 



from which we see that evaporating PBHs must form very 
early indeed in the history of the universe, in particular 
long before the epoch of nucleosynthesis (T ~ 1 MeV) at 
which the standard big bang evolution is well validated. 

In the 'standard' cosmology, the universe has been ra- 
diation dominated ever since the end of the reheating 
period after a phase of inflation at extremely high ener- 
gies, which was responsible for the generation of density 
perturbations. Under this assumption, the limits on the 
PBH density can be extrapolated backwards to the time 
of formation to give limits on the initial mass fraction 
of PBHs, fa = p p bha/ptot,i where p pb h,i and p tot<i are 
the PBH and total energy densities respectively, at the 
time at which the PBHs are formed. Substantial work 
has been done under this assumption However, 
there is no direct evidence requiring that the universe be 
radiation dominated at the high temperatures under con- 
sideration, and these limits can be greatly altered if the 
evolution of the Universe is more complex. In this paper 
we consider one of the most dramatic possible changes 
to the early evolution — the effect of a second period of 
inflation. 

A recent extension to the standard cosmology is ther- 
mal inflation Q|, a second period of inflation due to a 
scalar field known as a flaton. This has nothing to do 
with the normal period of inflation, which is still assumed 
to occur at some higher energy scale to solve the flatness 
and horizon problems and to generate density pertur- 
bations. Flaton fields, which are a consequence of su- 
persymmetric theories, have a vacuum expectation value 
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(vev) M ^> 10 3 GeV, even though their mass m is only of 
order the supersymmetry scale, ~ 10 2 to 10 3 GeV. Their 
potential is therefore almost flat. In the early universe 
these fields are held at zero by finite temperature effects, 
with false vacuum energy density Vb ~ m 2 M 2 . Once the 
temperature falls below V^ 4 , the false vacuum energy 
density dominates the thermal energy density. This false 
vacuum energy begins to drive a new period of inflation. 
This inflation continues until the temperature drops to 
T ~ Tn, at which point thermal effects are no longer 
strong enough to anchor the flaton in the false vacuum. 
The most popular possible identity for flaton fields are 
the moduli fields in superstring theory. 

A modulus field with a vev of order the Planck mass 
mpi (as is expected if the vev is non-zero) would produce 
many particles which would not decay before nucleosyn- 
thesis and hence would destroy the standard model of 
cosmology; this is the 'moduli problem' j^j. It cannot 
be solved in the same way as the monopole problem by 
invoking an early epoch of inflation lasting upwards of 60 
e-foldings, since the energy scale at the end of inflation 
is normally greater than 10 12 GeV and the moduli would 
be regenerated after inflation. To avoid too many moduli 
being regenerated requires § 

^ 4 <(10 7 tol0 8 )Gev(i^) 1/4 , (2) 

where Tr is the reheat temperature. A single period of in- 
flation at such a low energy scale would not be capable of 
producing the observed density perturbations, but a pe- 
riod of thermal inflation can solve this problem provided 
M is within a range of several orders of magnitude around 
10 12 GeV. Taking M ~ 10 12 GeV gives V 1/4 ~ 10 7 GeV 
so that around ln(10 7 /10 3 ) ~ 10 e-foldings of thermal 
inflation occur, sufficient to dilute the moduli existing 
before thermal inflation but small enough to not affect 
the density perturbations generated during the first pe- 
riod of inflation.^ 

The main effect of the period of thermal inflation is 
to dilute the density of PBHs, relative to radiation, by a 
factor of (af/ai) 3 ~ (10 4 ) 3 , where a; and at are the scale 
factors immediately before and after thermal inflation. 
In addition, there are two more subtle changes from the 
standard scenario. Firstly, for PBHs which form in the 
time between the two periods of inflation, the comov- 
ing scale to which a given mass of PBH corresponds is 
changed. During thermal inflation the horizon mass re- 
mains constant whilst the comoving Hubble radius grows, 



* Originally the flaton was taken to be the GUT Higgs [Q] 
so that M ~ Mgut ~ 10 16 GeV and V 1/4 ~ 10 10 GeV, 
leading to 15 e-foldings of inflation. However, successful nu- 
cleosynthesis requires M < 10 12 GeV and thermalization of a 
stable lightest supersymmetric particle, if one exists, requires 
M < 10 10 GeV §0. 



so a given PBH mass corresponds to a shorter comoving 
length scale and hence to a later stage in the original, 
density perturbation generating, epoch of inflation. Sec- 
ondly, thermal inflation introduces a missing mass range 
of PBHs. This corresponds to those comoving scales 
which enter the horizon before thermal inflation (pos- 
sibly forming black holes as they do so), and are then 
pulled back outside again during thermal inflation. Any 
new density perturbations might be expected to be small, 
since the energy scale of thermal inflation is much lower 
than the original inflationary period, and hence unable 
to form black holes when they re-enter the horizon again 
after thermal inflation.^] From Eq. (Ej), this corresponds 
to masses in the range 10 18 g < M ^10 26 g. 

Finally, PBH formation requires Sp/p on the relevant 
scales to be two or three orders of magnitude larger than 
the value required by COBE, which applies at very large 
scales. PBHs can therefore only be formed in significant 
numbers if the spectral index n of the density pertur- 
bations is significantly above unity, normally referred to 
as a blue spectrum. Limits on PBH formation allow up- 
per bounds to be placed on n; for the standard cosmology 
this was done by Carr et al. jl. We shall first re-examine 
their calculation, and then generalize it to include the 
possibility of thermal inflation. 

II. PRIMORDIAL BLACK HOLE FORMATION 

In order for a PBH be formed, a collapsing over- 
dense region must be large enough to overcome the 
pressure force resisting its collapse as it falls within its 
Schwarzschild radius. Consider a spherically symmet- 
ric region with density p greater than that of the back- 
ground, whose evolution will be governed by the positive 
curvature Friedmann equation uM. The perturbed re- 
gion stops expanding when H = 0, at which time the 

— 1/2 

region has size R c w S i i?; where i?i is its size at some 
arbitrary initial time and dj = (p~j — pi)/ pi is the initial 
density perturbation If the perturbed region con- 

tains enough matter to overcome any pressure forces, it 
will continue to contract. This requires that its radius 
exceeds the Jeans length, R c > i?j can s = c s t c , where c s 
is the sound speed. In a radiation-dominated universe, 

-Rjcans ~ tc/VS- Now R c /t c w RiftiS^ 2 , and since this 
expression is constant with time we can evaluate it at 
horizon crossing [R = t) leading to a constraint on the 
perturbations at horizon crossing 5 > 1/3. There is also 
an upper limit of 8 < 1; a perturbation which exceeded 
this value would correspond initially to a separate closed 



'Note however the standard perturbation calculation breaks 
down as the hypothesis that the initial state is the vacuum 
probably cannot be justified. It is not clear how to make the 
necessary generalization. 
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universe [|13j|5|, which is an inconsistent initial condition 
for our purposes. So for PBH formation we require the 
initial fluctuations to satisfy 



1/3 < 5< 1 



(3) 



When a perturbation satisfying the above condition 
crosses the horizon, a PBH will be formed with mass H 



M = 7 3/2 M H 



7' 



3/2 



1/2 
S'+form 



i 

fpi 



mpi , 



(4) 



where the background equation of state is p = jp, 7 being 
1/3 in a radiation-dominated universe. Here 



Af H = fp(i/- 1 ) 3 



(5) 



and g+ is the effective number of massless degrees of free- 
dom at this time. We define $ to be the initial mass 
fraction of PBHs, which is given by the fraction of the 
Universe satisfying Eq. (0): 



Ppbh,i 
Ptot,i 



1/3 



p(S)dS. 



(6) 



where p(8) is the probability distribution for S. 

Normally (for instance in large-scale structure studies) 
when one considers perturbations the probability distri- 
bution is assumed to be gaussian. This is well justified 
when the perturbations are small. However, Bullock and 
Primack |ll[] have recently challenged this assumption for 
PBH formation, since the perturbations cannot be very 
small if a significant formation rate is to be obtained. 
They typically find a suppression of large perturbations 
relative to the gaussian hypothesis. While this suppres- 
sion can be very dramatic when expressed in terms of 
the number of black holes formed, it actually does not 
lead to a large change in constraints on the perturbation 
spectrum. We shall therefore maintain the gaussian as- 
sumption for our derivations, and in Section [v| we shall 
assess the changes non-gaussianity introduces. 

In order to examine specific mass ranges, we have to 
smooth the density distribution, which is done in the nor- 
mal way using a window function W(kR), which we take 
to be a top-hat. For gaussian distributed fluctuations, 
the probability distribution of the smoothed density field, 
p(S(M)), is given by 



p(6(M))dS(M) 



1 



2trj(M) 



exp 



S 2 (M) 
2a 2 (M) 



d6(M) . 



Here a(M) is the mass variance evaluated at horizon 
crossing, defined for example in Ref. fl4| ] 



<j 2 {M) 



2tt 2 



P{k)W 2 {kR)k 2 dk. 



(8) 



where P(k) = (|<5fc| 2 ) is the power spectrum, the 5k be- 
ing the coefficients when 8{x) is Fourier expanded. The 
power spectrum is usually taken to have primordial form 
P(k) oc k n for simplicity; in general there is no reason 
to expect inflation models to produce power-law spectra 
over as wide as range of scales as we will need to consider, 
though there are models which do. 

The mass fraction of black holes is given from the above 

by 



/?(M) 



1/3 v2tt <j{M) 



exp 



2a 2 (M) 



dS . (9) 



Since the integrand is a rapidly falling function in the 
regime of interest, dropping by a factor exp (—0.5) every 
time (5 is increased by c, this integral can be approxi- 
mated by evaluating the integrand at <5 = 1/3 and mul- 
tiplying by cr(M), leading to 



j3(M) « a(M) exp - 



1 



18cr 2 (Af) 



(10) 



Strictly speaking, this is the mass fraction in black holes 
of mass greater than M, but in practice (3{M) is such 
a rapidly falling function that these can be taken to all 
have the same mass M. 

The final step is to relate the mass scales to comoving 
scales during inflation. In the notation of Ref. |fl4fl , the 
initial spectrum of perturbations is <5j?j(fc) cx fc"~Mvhere 



(^7) (?) ™: 



(11) 



The quantity 5n stays constant on scales above the Hub- 
ble radius and is a good estimate of the rms density con- 
trast at horizon entry |fl4|| . During radiation domination 
the comoving Hubble radius H~ 1 /a is proportional to 
T -1 , so that a given scale k crosses within the Hubble 
radius when 



k- 1 



H- 



cq 



(12) 



where subscript 'eq' refers to quantities evaluated at 
matter-radiation equality. Meanwhile the horizon mass 
Mh, Eq. (|J), varies as T~ 2 so that, in the absence of 
thermal inflation, we get 



( M H 



" /cqV^H.eq 



1/2 



(13) 



(7) Substituting this into the expression for 5#(k) we obtain 



^hor(-M) = er hor (M eq ) 



M 



(l-n)/4 



(14) 



where M oq is the horizon mass at matter-radiation equal- 
ity. We stress that this equation refers to the dispersion 
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Constraint 


Range 


Reason 


Qcvap < 0.04 


10 y g < M < 10 ia g 


Entropy per baryon at nucleosynthesis 


a C va P < 10- Me (M/m P i) 


M = 5 x 10 14 g 


7 rays from current explosions 


Qcvap < 6 x 10- 10 (M/m P1 ) 1/2 


10 9 g < M < 10 11 g 


nil production at nucleosynthesis 


acvap < 5 x W- 29 (M/m P1 ) 3r2 


10 1U g < M < 10 11 g 


Deuterium destruction 


aovap < 1 x 10- 5U (Af/mpi) Y/2 


10 11 g < M < 10 13 g 


Helium-4 spallation 



TABLE I. Limits on the mass fraction of PBHs at evaporation. 



at horizon crossing, not at constant timej^] 

The lightest black holes to form are those which en- 
ter the horizon immediately after inflation. For simplic- 
ity we shall assume prompt reheating, and Eq. (|l]) then 
gives the minimum mass. Carr et al. Q examine some 
consequences of delayed reheating. 



III. LIMITS ON THE PBH ABUNDANCE 
A. At evaporation 

The observational constraints on the mass fraction of 
black holes at evaporation, a cvap (M) = Ppbh/Vrad are 
well known , and are listed in Table |. 

To interpret these we need to relate the black hole mass 
to their lifetime. Carr || parametrizes the results of Page 
fl5|| , which were found numerically by considering the 
number of species which a black hole of given mass can 
emit at a significant rate, to give the following relation 
between PBH mass and lifetime 



9 x 10 



-27 



f(M) 



M 

Tg 



sec . 



(15) 



where f(M) depends on the number of particle species 
which can be emitted and is normalized to 1 for holes 
which emit only massless particles. Note that the bulk of 
the evaporation always takes place near the initial tem- 
perature, so in this expression one only needs / at the ini- 
tial mass and not as a time- varying quantity. Considering 
the number of spin states available for a PBH to evapo- 
rate into at the present day (g e s,o)', there are two polar- 
izations of photon plus three neutrino species which each 
have two spin states and give a contribution 7/8 times 



that of the photons, since they obey Fermi rather than 
Bose statistics. Therefore g cSfi — 2 + (3 x 2 x 7/8) = 7.25 
so that f(M) = <7eff,o/7-25 and 



1.2 x 10 4 



(iL 

ffeff V m Pl 



(16) 



with the value of g e g at the time of evaporation being 
taken. For M > 4 x 10 9 g the temperature at evapo- 
ration is sufficiently low (< 10~ 4 GeV) that g e s has its 
present day value of 7.25. Eq. (|l6| ) is often quoted with- 
out the factor of 1.2 x 10 4 ; this is certainly non- negligible 
when calculating PBH lifetime although other approxi- 
mations often made when limiting the initial PBH abun- 
dance (e.g. g c g ~ 1, 7 ~ 1 and M ~ Mh in Ref. ||) 
appear to largely cancel this factor. 



B. Present-day PBH density 

From Eq. @, PBHs of mass M > 5 x 10 14 g will 
not have evaporated by the present day, but their initial 
abundance can be constrained from the fact that their 
present mass density must not overdose the universe 



n 



pbh,0 







pbh,cq 



< l 



(17) 



C. Present-day relic density 



It has been argued |16| that PBHs may not evaporate 
completely, as originally assumed, but instead leave a 
relic with mass M rc i ~ mp\. If this is the case the present 
mass density of relics, which will remain from all PBHs 
with initial mass M < 5 x 10 14 g, similarly must not 
overdose the universe, leading to 



* This disagrees with Ref. in which a different scal- 
ing cr(M) oc M^ 1- ™^ 6 was used. This arises from assuming 
M oc k~ [i with no time-dependence — i.e. that the comoving 
mass density is conserved. This is true for matter domination 
but not for radiation domination where the comoving mass 
density decreases with time. Because our scaling is stronger 
with n, our final constraints on n are tighter than theirs. 
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-^form 
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Tfo 



< 1. 



(18) 
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o 



V^UIlo LI eilllL 


R an (to 
XYeillJifc! 


a-, < 3 X lCT lb <W' J s/M) 


10 9 s < AI < 10 la s 


Oi < 3 x 1(T 2 ' 


M ~ 5 x 10 14 g 


Qi < 3 x l(n 17 (10 9 g/AO 1/2 


10 9 g < M < 10 11 g 


Qi < 3 x 10" 22 (M/10 10 g) 1/2 


10 10 g < M < 10 11 g 


a; < 3 x 1(T 21 (A//10 9 g) 5/2 


10 11 g < M < 10 13 g 


Qi < 1 x 1(T 19 (Af/10 15 g) 1/2 


M > 10 15 g 


Qi < 0.1(M/10 15 g) 3/2 


M < 10 15 g 



TABLE II. Limits on the initial mass fraction of PBHs, 
without thermal inflation. We define Qi = /3i/(l — $) for 
compactness. The final constraint is the relic constraint which 
only applies if relics are assumed. 



IV. LIMITS ON INITIAL MASS FRACTION OF 
PBHS 

A. Standard evolution of the universe 

In constraining the initial mass function of the black 
holes, one needs to assume an entire history for the uni- 
verse from the time of formation to the present. Tradi- 
tionally, it has been assumed that the universe was radia- 
tion dominated up until the recent matter-dominated era, 
and then the limits of Section [II can easily be evolved 
backwards in time in order to constrain the initial mass 
fraction of PBHs. The energy density in radiation di- 
lutes as p r ad oc a~ 4 , whereas that in PBHs decreases 



more slowly, p pb h 



,-3 



( Ppbh 
V Prad 



a(M) e 



Therefore 



cvap 



1/2 



/cvap — „ , 

1 — Pi \ tf< 

Using Eqs. (f§) and @, and taking ~ 100 

a(M) cvap = 3.2-^- (— 
I- Pi \mpi 



(19) 



(20) 



The gravitational constraint can be evaluated simply: 

3 



PPBH.cq — Ppbh.i 



T, 



cq 



Pi 7T f^ v ™ m/ i ( T ( 



1 - fa 30 



form/Ti4 I ± oq 

9* 1 fori 



and 



Ptot, cq — 2— g° q T c 4 , 



"30"* oq: 
where gl° rm ~ 100 and g° q ~ 3,| so that 



(21) 



(22) 



^Note that g t is the effective number of degrees of freedom 
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FIG. 1. The tightest limits on the initial mass fraction 
of PBHs, Qj . The relic constraint is shown as a dotted line, 
emphasizing that it is not compulsory. 
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(23) 



Taking t cq = t n^ = e.S/i^Gyr x (4 x 10~ 5 li- 2 ) 3 / 2 and 
using Eq. (§) 



n 



li0 = 6 .7xlO«-^r^V 



/2 



Pbh, u , ■ v M , 

so that the constraint f2 pD h.o < 1 leads to 

1/2 

Vmpi 



(24) 



(25) 



The calculation for the relic limit can be carried out 
identically leading to 



J* < 1.5X10- 

I- Pi \m?\ 



3/2 



(26) 



with the extra factor of {M/m-pij from Eq. (|l8|). 

The various limits on the initial mass fraction of PBHs 
are displayed in Table [n] and illustrated in Fig. |l|. 



B. With a period of thermal inflation 

We model the period of thermal inflation by assum- 
ing that at T = 10 7 GeV the energy density in radiation 
splits into two components, with one degree of freedom 
becoming the inflaton and the remainder staying as ra- 
diation: 



as far as cosmology is concerned, evaluated at the photon 
temperature. It is lower than g e g quoted earlier since the 
cosmic neutrino background is at a lower temperature than 
the microwave background. 



5 



H 

a 




time 

FIG. 2. A schematic of the variation of the comoving Hub- 
ble radius (H -1 /a) with time for the standard evolution of 
the universe (solid line) and with thermal inflation (dashed 
line). Points A and B correspond to the end of the original 
period of inflation in the standard evolution and with ther- 
mal inflation respectively. Thermal inflation begins at point 
C and finishes at D, after which time the comoving Hubble 
radii must coincide. Between D and E the scales which are 
entering the Hubble radius are doing so for the second time 
so that no PBHs are formed in this region. The values of the 
comoving Hubble radius and the horizon mass at these points 
are displayed in Table [En]. We denote the current comov- 
ing Hubble radius and horizon mass as (i? _1 /a)o and A/ho 
respectively. 



Prad /5 rad + , 



(27) 



where p^ — 7r 2 (10 7 GeV) 4 /30 is the false energy density 
of the flaton field which drives thermal inflation, i.e. 



rr.4 

30 y 



30 vy 



1)T 4 



(28) 



Thermal inflation proper then commences once the radi- 
ation has redshifted sufficiently, at 



10 7 GeV 
(.9* I) 1 / 4 ' 



(29) 



when > /5 ra d and continues until the flaton field rolls to 
its true vacuum state at T = 10 3 GeV. We then assume 
that reheating is efficient so that the universe is reheated 
to T t j and the subsequent evolution has its standard form. 
Inefficient reheating, normally modelled as matter dom- 
ination, would make little difference. The duration of 
thermal inflation is negligible in terms of the evapora- 
tion time-scale, so its dominant effect on PBH evolution 
is to dilute their energy density relative to that of radi- 
ation. The radiation energy density is the same before 



Point 



A 



C 



D 



Comoving Hubble Radius 



T 



10 a GcV 



10 7 GcV 



ff- 1 \ T 

10 3 GcV 



Horizon Mass 



M H ,o 



M H,0(^)' 



TABLE III. Comoving Hubble Radii and horizon masses 
at points on Fig. 0. 



and after thermal inflation (assuming efficient reheating) 
whilst the energy density of PBHs is diluted by a factor 
Pi/ Pi — (af/ai) 3 ~ (10 4 ) 3 , where Oj and Of are the scale 
factors immediately before and after thermal inflation. 

Fig. H illustrates the variation of the comoving Hubble 
radius both with and without thermal inflation, with the 
Hubble radius and Afjj at important points being given 



in Table II] 



Provided that the PBHs do not come to dominate be- 
fore thermal inflation, the limits on the initial black hole 
mass fraction from the evaporation constraints are sim- 
ply weakened by a factor of 10 12 . The condition for the 
universe to be radiation dominated before thermal infla- 
tion commences is 



1-/3 



A 



T ti 



i - A V T ti 



Ti 



form 



< l 



Using Eq. (0) this requires 



T^W < 6 x 10_12 v 
l-A \"ipi 



/ M 



1/2 



(30) 



(31) 



If the PBHs come to dominate at T > 10 7 GeV then 
thermal inflation can only commence once the false en- 
ergy density of the flaton field dominates the energy den- 
sity on PBHs, > /?pbh- This delays the start of infla- 
tion, so that a smaller number of e-foldings of inflation 
occur; however, the dilution of the PBH energy density 
relative to that of radiation remains the same. Dur- 
ing the intermediate period p$ remains constant while 
Pphh oc a~ 3 , so pphh is rapidly reduced to below p,p and 
the resulting constraints from Hawking radiation on the 
initial mass fraction are only slightly tighter than if ther- 
mal inflation commences at T ti = 10 7 /3 GeV. However 
these limits are tighter than Eq. (pl|), so that in fact 
PBHs with M > 10 9 g cannot be produced with suffi- 
cient abundance that they come to dominate the universe 
before thermal inflation. 

In the case of lighter PBHs, which are only constrained 
by the present-day relic density, the condition for ra- 
diation domination before thermal inflation commences, 



G 



o 





R ATI ffp 


Qi <3x 1CT 4 (10 9 g/M) 


10 9 g < M < 10 i3 g 


tti <3x lCT lb 


2 x 10 14 g < M < 5 x 10 14 g 


on < 3 x 1(T 5 (10 9 g/A/) 1 / 2 


10 9 g < M < 10 11 g 


Qi < 3 x 1(T 10 (A//10 10 g) 1/2 


10 10 g < M < 10 11 g 


Qi <3 x l(T 9 (Af/10 9 g) 5/2 


10 n g < Af < 10 13 g 


Qi < 1 x 10 11 (A//10 15 g) 3/2 


M < 10 15 g 


a; < 1 x 10" 7 (A//10 15 g) 1/2 


10 15 g < M < 10 18 g 


Qi < 1 x 1(T 19 (M/10 15 g) 1/2 


Af > 10 26 g 



TABLE IV. Limits on initial mass fraction of PBHs if 
thermal inflation occurs 



Eq. (^Tj) , is more constraining than the requirement that 
f^rel.o < 1- These light PBHs can therefore come to 
dominate before thermal inflation and delay its start, as 
discussed above, although the resulting constraint from 
^rel,o < 1 is virtually the same as when the universe is 
radiation dominated at the start of thermal inflation. 

PBHs with mass M < 10 18 g are formed before ther- 
mal inflation and their energy densities will therefore be 
diluted by thermal inflation, so that the gravitational 
constraint for 10 15 g < M < 10 18 g and the relic con- 
straint will be weakened by a factor 10 12 . This allows 
PBHs with M < 10 9 g to be produced with initial abun- 
dance Pi close to one, although not arbitrarily so since 
sufficient thermal inflation to dilute the present day relic 
density must occur. During inflation Mh remains con- 
stant before increasing as T~ 2 again after inflation; how- 
ever, until the temperature falls to T = 10 3 GeV once 
more the scales that are entering the Hubble radius will 
be doing so for the second time having first entered at 
T > T t i (and possibly forming black holes) before being 
inflated away again. There are therefore no new density 
perturbations present to collapse into PBHs, so whilst 
T falls from 10 7 GeV to 10 3 GeV after thermal infla- 
tion no PBHs form leading to a 'missing' mass range 
10 18 g < M < 10 26 g. The gravitational constraints on 
PBHs with mass M > 10 26 g, which form after thermal 
inflation, are unchanged. 

The various limits on the initial mass fraction of PBHs 
are displayed in Table IV and illustrated in Fig. |[ 




Excluded 
mass range 



10 



15 



30 



25 



log M 

FIG. 3. The tightest limits on the initial mass fraction 
of PBHs Qi if thermal inflation occurs, on the same vertical 
scale as Fig. |l|. The gap 10 ls g < M < 10 26 g is the excluded 
mass range. For M < 10 9 g a large initial mass fraction ~ 1 
of PBHs is allowed. 



When inserted in Eq. (^) and numerically integrated, 
this allows us to normalize a(M) for a chosen fixed M. 
We take the COBE data to correspond to a scale equal 
to the present Hubble radius i?o = 3000/i _1 , so that 
Mo = 1 x 10 56 g is the present horizon mass which gives 
a normalization cr(10 56 g) = 9.5 x 1CT 5 . To the level of 
accuracy at which we are working <r(10 56 g) varies only 
slowly with n. Putting this in Eq. ([L4|) gives 



cr hor (M) = 9.5 x 10" 




i)/4 



(33) 



Here we ignore the change in slope to (1 — n)/6 which 
occurs at matter domination (corresponding to a mass 



Mo t oq /t 



10 



which only changes the constraint 



on n by about 0.01. 

Taking logarithms in Eq. ( |l0| ) gives 

a(M) = 0.15 (log 10 a(M) - log 10 ft) 



-1/2 



(34) 



Since PJ(l - p { ) <C 1, we can take P i w P-J{\ - Pi) and 
find the upper bounds on n corresponding to each of the 
limits found in the previous section. 



V. LIMITS ON THE SPECTRAL INDEX 



Using Eq. (14), the limits on the initial mass fraction 
can be used to constrain the spectral index of the den- 
sity perturbations n. The four-year fitting function to the 
COBE data, assuming negligible contribution from grav- 
itational waves, gives the normalization at the present 
Hubble scale (k = a H ) [p7l 



<5 H (n) = 1.91 x 10~ b exp [1.01(1 - n)] 



(32) 



A. Standard evolution of the Universe 

The appropriate constraint is normally determined by 
the lightest PBHs that can form, given by Eq. (|l|). The 
tightest limit is n < 1.22 from the deuterium destruc- 
tion constraint evaluated at M ~ 10 11 g, although all 
the constraints due to the evaporation of the PBHs re- 
quire n < 1.24. The tightest constraint from the limit on 
the present density of PBHs is n < 1.31 at M ~ 10 15 g, 
weakening with increasing M . The relic constraint may 
place a stronger constraint on n, depending on the min- 
imum mass of PBHs produced, which is determined by 
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1 RH ( GeV * 

FIG. 4. The variation of the limits on n with reheating 
temperature from the relic constraint (lower line) and from 
<5(Mmm) < 1. 



the reheating temperature Tpjj after the initial period of 
conventional inflation 



M„ 



Mn 



RH 



(35) 



If Trh < 10 9 GeV, there is not sufficient time for ppbh 
(and after evaporation p re i) to increase by enough rela- 
tive to /5 ra d for the relics to dominate today. For com- 
parison, in Fig. ^| we plot the variation of the limit on 
n with reheating temperature and the limit from the 
simple requirement that (5(M m ; n ) < 1, where M m - U1 = 
(Tpi/Trh) 2 tIpi is the mass of the lightest PBHs formed 
immediately after the first period of inflation. 

In summary, the tightest Hawking radiation constraint 
on n we obtain is n < 1.22, although if the reheating tem- 
perature is sufficiently high the existence of relics may 
lead to a tighter limit. In Ref. Q a much weaker limit of 
n < 1.48 from Hawking radiation constraints was found, 
which becomes tighter if relics are formed (n < 1.4 if 
Trh = 10 16 GeV). Our stronger constraint arises from 
our correcting of two errors in their paper, which both 
go the same way. The first is that we use the correct scal- 
ing law a{M) oc M^ 1_n ^ 4 for the variance at the horizon 
scale during radiation-domination, as mentioned earlier. 
The second is that our normalization to COBE is much 
higher. They omitted a numerical prefactor -\/5127r/75 
in the power spectrum expression, and also assumed that 
the normalization of 5u(k) (which from the four-year 
COBE data is 2 x 10 -5 ) and a(M) were interchangeable. 
In combination this raises the COBE normalization by 
a factor of over twenty. A much smaller additional cor- 
rection is that the normalization from the COBE four- 
year data |l^] is higher than that from the first year's 
data Il9). 



B. With a period of thermal inflation 

During thermal inflation the comoving Hubble radius 
(H^ 1 /a) varies as T so Eq. (|l2|), for when a given co- 
moving scale k crosses the Hubble radius, becomes 



cq 



whilst the relation between horizon mass and tempera- 
ture remains unchanged, since during inflation the hori- 
zon mass is constant so that for PBHs formed in between 
the two periods of inflation a given mass PBH will cor- 
respond to a larger scale than in the standard scenario. 
This leads to the modification of Eq. ( [14} ) : 



(Thor(M) = O-hor(Meq) 



M /l0 3 GeVV 



M, 



eq 



V 



(l-n)/4 



(37) 



In this case the tightest constraint is n < 1.29 from 
the deuterium destruction constraint evaluated at M ~ 
10 10 g, with all the constraints due to evaporation requir- 
ing n < 1.34. In this case the constraints from gravitation 
and relics are weaker. The limit from the present-day 
density leads ton< 1.35 for M ~ 10 15 g PBHs, which 
are diluted by the thermal inflation, and n < 1.49 for 
M ~ 10 26 g PBHs, which are the lightest formed after 
thermal inflation. For Trh < 10 14 GeV the relics do not 
constrain n, since even if /3; is close to one the relics will 
be diluted away. However, for Trh = 10 16 GeV we find 
n < 1.3 is required although this limit rapidly weakens 
as the reheat temperature falls towards 10 14 GeV. 



C. The effect of non-gaussianity 

We now return to the issue of the gaussianity assump- 
tion used to obtain Eq. (Q). Bullock and Primack |Tl| ] 
have stressed that the normal justification of gaussian- 
ity relies on the perturbations being very small, some- 
thing which can no longer be justified when considering 
PBH formation. Unfortunately, the non-gaussian cor- 
rection is strongly model dependent, and in detail must 
be examined case-by-case. They numerically study three 
different 'toy' models, in one case finding negligible non- 
gaussianity but in the other two finding a very significant 
suppression in the number of large perturbations, which 
are of course exactly those utilized for PBH formation. 

Expressed in terms of the probability of high density 
perturbations, the suppression can be very dramatic; in 
one of their toy models the 6-sigma perturbations are 
suppressed by a factor of 10 150 ! However, despite that 
the effect on the constraint on n is not large, because 
that asks a rather different question, namely how much 
larger does the variance a(M) have to be so that the 
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non-gaussian perturbations reproduce the number den- 
sity of gaussian ones? In their most extreme example, 
the answer is about three times; the perturbations with 
the appropriate number density correspond to about 3- 
sigma perturbations in the non-gaussian case rather than 
the 9-sigma or so perturbations of the gaussian case [jllj . 

The conclusion then is that non-gaussian effects are 
model-dependent, and in the worst tested case weaken 
the constraints on a\ ioI (M) by about a factor three. Non- 
gaussianity clearly cannot do much more than this, as the 
low required number density keeps us to the tail regard- 
less of the amount of non-gaussianity. From Eq. (14), 
using the COBE normalization to keep crhor(-Mo) fixed, 
this weakens the constraint on n, in the worst case of 
non-gaussianity, by about 0.05. 



VI. CONCLUSIONS 

The interpretation of constraints on primordial black 
holes depends sensitively on the whole history of the uni- 
verse from their time of formation. We re-examined the 
constraints assuming the standard radiation-dominated 
cosmology, and by correcting two errors in Carr et al. || 
found a significantly tighter constraint on the spectral 
index than they did, n < 1.25. This is in fact presently 
the tightest constraint on the spectral index, being some- 
what stronger than large-scale structure constraints in 
the most general cosmologies, and much stronger than 
the constraint from distortions to the microwave back- 
ground spectrum po[ . On the other hand, its application 
requires a belief that the spectral index remains constant 
over a much wider range of scales than the others, which 
is certainly possible but not mandatory. In general cir- 
cumstances one must impose the general constraints on 
the formation rate we derived in Section IV. 



We have analyzed the changes to the standard sce- 
nario brought about if a period of thermal inflation takes 
place in the early universe after the black holes form. 
Thermal inflation leads to a significant weakening on 
the constraint on the density perturbation spectrum. In 
Section IV we recomputed the constraints on the initial 



mass fraction of black holes; especially at low masses 
the constraints become very weak indeed. In Section 
we assumed a power-law spectrum of perturbations and 
constrained the spectral index n; we found the con- 
straint weakened (relative to the standard cosmology) 
to n < 1.3. A novel additional feature is that thermal 
inflation predicts a missing mass range for black holes, 
extending up from 10 18 g to 10 26 g. It will be hard to 
probe this range as black holes of these masses have neg- 
ligible evaporation. However, if for some reason thermal 
inflation can start at a higher energy than currently sup- 
posed, say 10 10 GeV, then the missing mass range could 
extend down into the evaporating regime. 

It seems quite likely that density perturbations large 
enough to form PBHs will exhibit a significantly non- 



gaussian probability distribution, as emphasized by Bul- 
lock and Primack jll]]. However, we have shown that this 
does not much alter the constraints; in the worst case it 
corresponds to a weakening of about 0.05 (regardless of 
whether or not thermal inflation occurred) which is com- 
parable to the change due to thermal inflation. But typ- 
ically the correction is expected to be quite a bit smaller 
than this |p"l|| , and hence smaller than the uncertainty in 
the cosmological model at early times. 
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